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1. Introduction 

For spherically symmetric accretion onto a non-rotating compact object, a novel 
mechanism was proposed by Protheros and Kazanas,^ (PK83 hereafter) and Kazanas 
and Ellison"^ (KE86 hereafter) , in which the kinetic energy of the accreting mate- 
rial was shown to be randomized by incorporating a steady, collision-less, rela- 
tivistic hadronic-pressure-supported spherical shock surface around the accreting 
Schwarzschild black hole which produces a non-thermal spectrum of relativistic 
protons. Recently Das'* (D99 hereafter) has explicitly calculated the exact location 
(radial distance measured from the central accretor in units of Schwarzschild radius 
fg — ^^^^j^, where AIbh is the mass of the black hole, G is the Universal gravi- 
tational constant and c is the velocity of light in vacuum) of the above mentioned 
shock in terms of only three accretion parameters, namely, the specific energy of the 
flow £, accretion rate Msdd (scaled in units of Eddington rate) and the adiabatic 
index jin of the inflow for spherically-symmetric, transonic accretion of adiabatic 
fluid onto a Schwarzschild black hole. By solving the set of hydrodynamic equations 
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describing the motion of accreting material under the influence of modified Newto- 
nian potential proposed by Paczyiiski and Wiita,^ it has been shown there in D99 
that it is possible to construct a self-consistent inflow-outflow system where a part 
of the accreting material may be blown as wind from the spherical shock surface 

and the mass outflow rate Rm (the measure of the fraction of accreting material 
being 'kicked off' as wind) was computed in terms of various accretion as well as 
shock parameters. 

However, other than the Paczynski-Wiita potential, a number of modified New- 
tonian potentials of various forms are also available in the literature which accu- 
rately approximate some general relativistic features of rotating accretion around 
Schwarzschild black holes. Such potentials may be called 'pseudo - Schwarzschild' 
potentials because they mimic the space-time around a non-rotating/slowly rotating 
compact object. Recently Das and Sarkar*^ (DS hereafter) has shown that though 
these so called pseudo potentials were originally proposed to mimic the relativistic 
effects manifested in disc accretion, it is quite reasonable to use most of these poten- 
tials in studying various dynamical and thermodynamic quantities also for spherical 
accretion on to Schwarzschild black holes and it was established that along with the 
Paczynski and Wiita potential,^ two other potentials proposed by Artemova et. al.^ 
also provide reasonably good approximation to the complete general relativistic de- 
scription of transonic, spherically symmetric accretion on to a Schwarzschild black 



Remembering that the free-fall acceleration plays a very crucial role in Newtonian 
gravity, Artemova et. al.^ proposed the following two potentials to study disc ac- 
cretion around a non-rotating black hole. The first potential proposed by them 
produces exactly the same value of the free-fall acceleration of a test particle at 
a given value of r as is obtained for a test particle at rest with respect to the 
Schwarzschild reference frame, and is given by 



The second one gives the value of the free fall acceleration that is equal to the value 
of the covariant component of the three dimensional free-fall acceleration vector of 
a test particle that is at rest in the Schwarzschild reference frame and is given by 



For Eq. 1 (a-b) , r represents the radial co-ordinate scaled in the unit of Vg . Efficien- 
cies produced by $i and $2 are —0.081 and —0.078 respectively. As both $1 and 
$2 stems from the consideration of free fall acceleration and calculates the depen- 
dence of free fall acceleration on radial distance in the Schwarzschild metric (which 
describes a spherically symmetric gravitational field in vacuum), it appears to be 
quite justified to use those potentials to study spherically symmetric accretion. 
Owing to the fact that $1 and $2 may be used to mimic the spacetime around a 



hole. 




(la) 




(16) 
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spherically accreting non-rotating black hole quite nicely, we believe that along with 
the calculation presented in D99, it is equally important to investigate the various 
features of the accretion powered spherical outflow using these two potentials. In 
this paper we precisely do this, first we formulate and solve the equations governing 

the inflow and outflow using $i and $2 (in presence of a steady, standing spherical 
shock as described in previous paragraphs) and then self-consistently connect the 
accretion and wind solutions to calculate what fraction of the accreting material is 
being blown as wind. Also we study the dependence of this fraction on various flow 
parameters. The plan of this paper is as follows. In next section we describe how 
to formulate and solve the governing equations. In §3, we present our results. In 
§4, the results will be reviewed and conclusions will be drawn. 



2. Governing equations and solution procedure 

2.1. Inflow model 

Hereafter we will refer to these two potentials as <I>i in general where {i = 1,2} 
would correspond to $1 (Eq. (la)) and $2 (Eq- (lb)) respectively. We assume 
that a Schwarzschild type black hole spherically accretes fluid obeying a polytropic 
equation of state. The density of the fluid is p(r), r being the radial distance 
measured in the unit of Schwarzschild radius Vg . We also assume that the accretion 
rate (scaled in the unit of Eddington rate Mecu) is not a function of r and we 
ignore the self-gravity of the flow. For simplicity of calculation, we choose the 
geometric unit where the unit of length is scaled in units of r^, units of velocity in 
units of c . All other physically relevant quantities can be expressed likewise. We 
also set G = c = 1 in the system of units used here. It is to be mentioned here 
that one fundamental criterion for formation of hydrodynamic outflows is that the 
outflowing wind should have a positive Bernoulli constant which means that the 
matter in the post-shock region is able to escape to infinity. However, positiveness 
in Bernoulli's constant may lead to another situation as well where shock may 
quasi-periodic-ally originate at some certain radius and propagate outwards without 
formation of outflows. So formation of outflows is one of the possible scenarios 
when we focus on the positive energy solutions. In this paper we concentrate only 
on solutions producing outflows. Another assumption made in this paper is to 
treat the accreting as well as post-shock matter as a single temperature fluid, the 
temperature of which is basically characterized by proton temperature. 
For any pseudo-Schwarzschild potential the equation of motion for spherically 
accreting material onto the accretor is given by 

du du 1 dP ^ ' „ ,^ . 

+ + =0 2a 

at or p or 



where symbols have their usual meaning. The first term of Eq. 2(a) is the Eulerian 
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time derivative of the dynamical velocity at a given r, the second term is the 'advec- 
tive' term, the third term is the momentum deposition due to pressure gradient and 
the final term is due to the gravitational acceleration for a particular ith potential 
The continuity equation can be written as 



dp , 1 d ^ 2\ 







{2b) 



For a polytropic equation of state, i.e., p = Kp^ , the steady state solutions of eqn. 
2(a) and eqn. 2(b) are 

1) Conservation of specific energy £ of the flow: 



2 7m - 1 



where 7i„ is the adiabatic index of the inflow (accretion) and 
2) Conservation of Baryon number (or accretion rate Mi„): 



Mi, 



AiTpur 



Using Min as the entropy accretion rate where M. = Min'Jin'^*" 
3(b) can be rewritten as (see DS and references therein): 



Mi 



Aira-^in-^ur 



(3a) 

(36) 
eqn. 

(3c) 



It is now quite straightforward to derive the spatial gradient of dynamical velocity 
(^) . for flow in any particular ith potential $i as 



^i 



u - ^ 



(4a) 



where 



denotes the absolute value of the space derivative of $i, i.e.. 



dr 



Since the flow is assumed to be smooth everywhere, if the denominator of Eq. 
4(a) vanishes at any radial distance r, the numerator must also vanish there to 
maintain the continuity of the flow. One therefore arrives at the so called 'sonic 
point conditions' by simultaneously making numerator and denominator of Eq. 4(a) 
equal to zero and the sonic point conditions can be expressed as follows 



where superscript i denotes the specific value of sonic quantities for a particular ith 



potential and 



is the value of (^) evaluated at the corresponding sonic 
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point r^. The value of sonic point rj. for any ith potential $i can be obtained by 
algebraically solving the following equation 



£- 



1 / 7m + 1 



1 



2 



= 



(4c) 



where 



is the value of ith potential at the corresponding sonic point r* . Similarly, 



the value of (^). for any i>j at its corresponding sonic point rl can be obtained by 
solving the following quadratic equation: 



+ 2.829{-n„-l) 



+ {'^lin - 1) 



= 



{Ad) 



where 



is the value of ^p^- at the corresponding critical point r]. 



One can simultaneously solve eqn. 3(a) and eqn. 3(b) (alternatively, eqn. 3(a) and 
eqn. 3(c)) for any specific for a fixed value of £ and 7i„ to obtain various dynam- 
ical and thermodynamic flow quantities. In this work we normally use the value of 
7j„ to be I . Though far away from the black hole, optically thin accreting plasma 
may not be treated as ultra-relativistic fluid (by the term 'ultra-relativistic' and 
'purely non-relativistic' we mean a flow with 7 = | and 7 = § respectively, accord- 
ing to the terminology used in Frank et. alJ) in general, close to the hole it always 
advects with enormously large radial velocity and could be well-approximated as 
ultra-rclativistic flow. As because our main region of interest, the shock formation 
zone, normally lies close to the black hole (a few tens of Vg away from the hole or 
sometimes even less, see results and figures in §3), we believe that it is fairly justifi- 
able to assign the value | for 7j„ in our work. However, to rigorously model a real 
flow without any assumption, a variable polytropic index having proper functional 
dependence on radial distance (i.e., 7i„ = ^mir)) may be considered instead of 
using a constant 7i„, and equations of motion may be formulated accordingly, which 
we did not attempt here for the sake of simplicity. Nevertheless, we keep our option 
open for values of 7j„ other than | as well and investigated the outflow for an range 
of values of 7i„ for a specific value of £ and Msdd (Fig. 4, §3.3). The same kind of 
investigations could be performed for a variety of values of £ and MEdd and a set 
of results may be obtained with a wide range of values of 7i„ which tells that our 
calculation is not restricted to the value 7j„ = | only; rather the model is general 
enough to deal with all possible value of 7i„ for polj^ropic accretion. 
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2.2. Shock parameters and the outflow model 

As already mentioned, a steady, collision-less shock forms due to the instabilities 
in the plasma flow (see D99 and references therein for details of shock formation 
mechanism). We take this shock surface to be the effective physical barrier around 
the black hole which might be responsible to generate the outflow. We assume 
that the presence of a collision-less steady standing spherical shock discussed in 
this work may randomize the directed in-fall motion and at the shock surface the 
individual components of the total energy of the flow (which is a combination of the 
kinetic, thermal and gravitational energy) gets rearranged in such a manner that 
the thermal energy of the post-shock matter dominates (due to enormous shock 
generated post-shock proton temperature) over the gravitational attraction of the 
accretor and a part of the in-falling material is driven by thermal pressure to es- 
cape to infinity as wind. In §3. we show that for any shock solution, the mass-loss 
rate normally co-relates with the post-shock proton temperature which essentially 
supports the validity of our assumption. We also assume the effective thickness of 
the shock Ash to be small enough compared to the shock standoff distance Vgh, and 
that the relativistic particles encounter a full shock compression ratio while crossing 
the shock. 

At the shock, density of matter will shoot up and inflow velocity will drop abruptly. 

If {p- , M_) and (p+, u+) arc the pre- and post-shock densities and velocities respec- 
tively at the shock surface, then 

^ - i? - — (5) 

P- u+ 

where Rcomp is the shock compression ratio. For high shock Mach number solution, 
the expression for Rcomp 

can be well approximated as 

Rcomp = lA4:Msh~* (6) 

where Mgh is the shock Mach number and Eq. (5) holds for Mgh ^ 4.0.^ 

In terms of various accretion parameters, shock location can be computed as (D99): 

where app is the collision cross section for relativistic protons, Ush and Mgh are the 
dynamical flow velocity and the Mach number attained at the shock location, Msdd 
is the mass accretion rate scaled in units of Eddington rate. One can understand 
that 

(ush, Msh) = C MEdd, 7m) (7a) 

where C, has some complicated non-linear functional form which cannot be evalu- 
ated analytically, but the value of Ush and Msh can easily be obtained in terms of 
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^£,MEdd,7in^ by numerically solving Eq. 3(a-b) and Eq. 7 (with the help of Eq. 
4(a-d)) simultaneously. Hence one can write 

rsh = (s, Msdd, 7m) (7&) 
where ^ has some functional form other than that of (. 

In ordinary stellar mass- loss computations (Tarafder,^ and references therein), the 
outflow is assumed to be isothermal till the sonic point. This assumption is probably 
justified, since copious photons from the stellar atmosphere deposit momenta on 
the slowly outgoing and expanding outflow and possibly make the flow close to 
isothermal. This need not be the case for outflow from black hole candidates. Our 
effective boundary layer, being close to the black hole, are very hot and most of 
the photons emitted may be swallowed by the black hole itself instead of coming 
out of the region and depositing momentum onto the outflow. Thus, the outflow 
could be cooler than the isothermal flow in our case. We choose polytropic outflow 
with a different polytropic index 7^ < 7j„ due to momentum deposition. In our 
calculation we also assume that essentially the post-shock fluid pressure and the 
post-shock proton temperature controls the wind formation as well as the barionic 
matter content of the outflow. 

The adiabatic post-shock sound speed a^^ and the post-shock temperature T^"^ 
(which is basically the temperature of the protons according to our one-temperature 
fluid approximation) can be calculated as: 



at, = (8a) 

V Psh 

and 

TX = ^ (86) 

where and pj^ are the post-shock pressure and density of the flow at shock 
location Vsh respectively. For low energy accretion ('cold' inflow, so to say) which 
is appropriate to produce a high shock Mach number solution, one can assume that 
the pre-shock thermal pressure (p^ft.) ™-^y be neglected compared to its post-shock 
value (pth) and to the pre-shock ram pressure (u^^)^). One can obtain the 

value of using Eq. (5-6) and from the total pressure balance condition at shock 
as, 

= «.)^«. (%^) (9) 

Combining Eq. (5-9), post-shock sound velocity and temperature obtained at the 
shock surface can be rewritten as: 



7o /1.44MA-1\ , , 



\ ^Edd \ m\ 
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and 



M, 



(106) 



sh 



General form of the conservation equations governing the polytropic outflow will be 

the same as Eq. (3a - 3c) with different polytropic index and total specific energy 
and entropy (f > £ , M.in < M and 7i„ > 7o, see D99). So we can write: 



= + + 

2 7o - 1 

Mont = 47rp°M°r" 



and 



(11a) 
(116) 

(11c) 

where where £ is the specific energy of the outflow which is also assumed to be 
constant throughout the flow and Jvi = Moutlo '<° ^ ^ K° - 1 is the entropy accretion 
rate of the outflow. 70 < 7ira as already mentioned. Any sub/ super-script indicates 
that the quantities are measured for the outflow. Like Eqs. (4a - 4d), one can easily 
write the sonic point conditions and the velocity gradient of the outflow as: 



M =47r(a°)~uV 



du° 
dr 



u" - 



1 flo + l 



\ 



(rM° , 



and 



(12a) 

(126) 
(12c) 



rlii° 

(l + 7o)(^) +2.829(7o-l)' 



+ (27o - 1) 



We then define the mass outflow rate Rm as: 



Mir,. 



du° 



=0 



(12d) 



(13) 
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It is obvious from the above discussion that Rm should have some comphcated 
non-hnear functional dependence on the following accretion and shock parameters: 

Rm = '^ Msdd, Tsh, Msh, Rcomp, 7m, 7o) (14a) 

As rsh,Msh and Rcomp can be found in terms of £,MEdd and only, ultimately 
it turns out that: 

Rm = ^ (S, MEdd, lin, 7o) (14&) 

Where ft has some complicated functional form which cannot be evaluated analyt- 
ically. 

2.3. Simultaneous solution of inflow-outflow equations 

In this work, wc arc interested in finding the ratio of Mout to Mm (Eq- (13)), and 
not the explicit value of Mout- Also note that the primary goal of our present work 
was to compute the outflow rate and to investigate its dependence on various inflow 
parameters but not to study the coUimation procedure of the outflow. 
Before wc proceed in detail, a general understanding of the transonic inflow out- 
flow system in the present case is essential to understand the basic scheme of the 
solution procedure. Let us consider the transonic accretion first. Infalling matter 
becomes supersonic after crossing a saddle- type sonic point, the location of which is 
determined by Msdd, 7mi|- This supersonic flow then encounters a shock (if 
present), location of which {rgh) is determined from Eq. (7). At the shock surface, 
part of the incoming matter, having a higher entropy density (because shock in a 
fluid flow generates entropy), is likely to return as wind through a sonic point other 
than the point through which it just entered. Thus a combination of transonic 
topologies, one for the inflow and other for the outflow (passing through a different 
sonic point and following topology completely diff'erent that of the 'self-wind' of 
the accretion), is required to obtain a full solution. So it turns out that finding 
a complete set of self-consistent inflow outflow solutions ultimately boils down to 
locate the sonic point of the polytropic outflow and the mass flux through it. Thus 
a supply of parameters £,MEdd, lin and 70 make a self-consistent computation 
of Rm possible. Here 70 is supplied as free parameter because the self-consistent 
computation of 70 directly using £, MEdd and has not been attempted in this 
work; instead we put a constrain that 70 < 7m always and for any value of 7,^. In 
reality, 70 is directly related to the heating and cooling processes taking place in 
the outflow. 

We obtain the inflow sonic point Vc by solving Eq. (4c). Using the fourth order 



"('•) 

a(r) 



are com- 



the Runge Kutta method, u{r), a(r) and the inflow Mach number 
puted along the inflow from the inflow sonic point Vc till the position where the 
shock forms. The shock location is calculated by solving Eq. (7). Various shock 
parameters (i.e., density, pressure etc at the shock surface) are then computed self- 
consistently. 
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For outflow, with the known value of £' and jo, it is easy to compute the location of 
the outflow sonic point r° from Eq. (12c). At the outflow sonic point, the outflow 
velocity u° and polytropic sound velocity o° is computed from Eq. (12b). Using 
Eq. (12a) and (12d), (^) and (^)^ is computed as was done for the inflow. 
Rungc -Kutta method is then employed to integrate from the outflow sonic point 
r° towards the black hole to find out the outflow velocity u° and density p° at the 
shock location. The mass outflow rate Rj^ is then computed using Eq.(13). 

3. Results 

3.1. Shock location as a function of £ and Msdd and related post-shock 
quantities 

For a particular value of £, Msdd and 7i„, the shock location (measured from the 
black hole in units of r^) can be calculated using Eq. (7). As Ugh and Mgh is 
a function of £ , Msdd, and 7i„, Vsh will also change with the change of any of 
these accretion parameters. In figure 1, we show the variation of Vsh as a func- 
tion of £ for three different values of Msdd- While £ is plotted along the X axis, 
shock location (in logarithmic scale) is plotted along the Y axis for a fixed value of 
7m (= |)- Three different curves drawn by solid lines for $2 and dotted lines for 
$1, are plotted for three different values of Msdd (=0.25, 1.0, 1.75). The lowermost 
curves (for each of the pseudo-potentials) corresponds to the value Msdd = 0.25 as 
shown in the figure. For both $1 and $2, other two curves from bottom to top, 
correspond to MEdd equal to 1.0 and 1.75 respectively. For any pseudo-potential 
one should note that different curves terminate (in the direction of increasing £) at 
different points which indicates that shock formation is not a generic phenomena, 
i.e., shock does not form for any value of £, MEdd, and 7i„, rather a specific region 
of parameter space spanned by M^dd, 7m| allows shock formation. Both sub- 
as well as super-Eddington accretion allows shock formation as shown in the figure. 
For any pseudo-potential, while the shock location non-linearly anticorrelates with 
£ (for a fixed value of MEdd and "fin), it correlates (non-linearly) with Msdd (for a 
fixed value of £ and 7i„). The maximum value of £ for which shock may form for 
a fixed value of Msdd and 7i„ , increases with increase of the accretion rate of the 
flow. It is observed that (figuratively not shown in the paper) the shock location 
also non-linearly anti-correlates with 7^ which means that for both and $2, 
non-relativistic super-Eddington accretion with low specific energy of the flow is 
a proper combination to produce the shock closest to the black hole (this result 
has significant importance in studying the amount of barionic load in the wind, see 
§3.3). Once the shock location is known, one can easily calculate any post-shock or 
shock related (the shock compression ratio for example) quantity using the corre- 
sponding equations derived in §2.2. One important physical quantity of our interest 
is the post-shock proton temperature of the flow (which is practically the temper- 
ature of the outflow according to our one-fluid approximation, see §2.1) which 
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can be computed using Eq. 10(b). We have seen that for any pseudo-potential, T^"^ 
correlates with £ as well as with 7j„, hence high energy pure-nonrelativistic flow 
would produce hotter outflow as well as it is evident that the closer the shock forms 
to the accretor, the higher becomes the post-shock proton temperature of the flow. 
Whatever observations arc presented in the above paragraph, is commonly applica- 
ble for both $1 and $2- However, there are a number of differences in magnitude of 
post-shock quantities (even in the shock location) observed when the flow is studied 
for two different pseudo-potentials. From figure 1, one can observe that for the 
same values of £, Msdd and 7j„, the shock forms relatively closer to the black hole 



for $2 compared to the case for $i. If rgh 



is the shock location obtained using 



any particular ith pseudo-potential, we see that: 



(15) 



for a fixed value of |f , MEdd, lin | • This deviation is more prominent for flows with 
higher specific energy and the following quantity 



decreases in the direction of low £. Also we note that the maximum value of £ for 
which the shock forms for a fixed value of Msdd and 7i„ (let us define that energy 



clS ^Tf 



for any ith pseudo-potential), is higher for flows in $2, i-e.. 



> £r, 



but unlike Svsh, S£max does not show any specific dependence on inflow parameters. 
It is also observed that (not shown in the figure) for a fixed values of energy and 
accretion rate, if we study the shock location as a function of the polytropic index 
of the infiow, the inequality presented in equation (15) is still maintained and Srsh 
decreases as the flow tends to its pure non-relativistic limit. As the post-shock 
proton temperature is inversely proportional with the shock location, we find 
that, 



for both the cases when: 

a) £ is being varied keeping Msdd and 7,, 



-^sh 



constant. 



(16) 
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b) Flow is being studied as a function of 7i„ for fixed values of S and Msdd- 
However, ST^ decreases with decrease in £ but with increase in 7i„. 

3.2. Combined integral curves of motion 

Figure 2a and 2b show two typical solutions which combines the accretion and 
outflow for flows in $i and $2 respectively. The accretion parameters used are £ = 

0.0003, Msdd =0.5 and 7^^ = | corresponding to ultra-rclativistic inflow. For both 
of the figures 2a and 2b, the solid curve AB represents the pre-shock region of the 
inflow and the solid vertical line BC with double arrow at rgh represents the shock 
transition. Shock locations (9.9 Vg for the Fig. 2a and 10.6 rg for the Fig. 2b) is 
obtained using the Eq. (7) for a particular sot of inflow parameters mentioned above. 
Three dotted curves show the three different outflow branches corresponding to the 
three different adiabatic indices 70 of the outflow. Prom left to right, the values of 
7o are 1.3, 1.25 and 1.2 respectively with respective mass-outflow rates as 1.558 x 
10-^ 7.1487 X 10-5 and 6.249 x 10"^ (for Fig. 2a) and 1.588 x lO'^, 7.47 x 10"^ 
and 6.765 x 10~^ (for Fig. 2b) respectively which indicates that for a given value 
of £, Msdd and 7^, i?,,, c;orrelates with 70. It is evident from the flgure that the 
outflow moves along the solution curves in a completely different way to that of 
the 'self- wind' solution of the inflow (solid line marked by CD in Fig. 2a and 2b). 
Also, the sonic points for all the outflowing branches are different to those of the 
accretion 'self- wind' system which is designated as Pg. While Pg = 1250.833 rg for 
for Fig. 2a and 1251.665 rg for Fig. 2b, the sonic points of the outflowing branches 
corresponding to 70 = 1.3, 1.25 and 1.2 are 1528.55 rg, 2084.033 r^, 2917.31 rg for 
Fig. 2a and 1529.32 rg, 2084.733 fg, 2917.95 rg for Fig. 2b respectively which 
indicates that the outflow sonic point increases with a decrease in the adiabatic 
index of the outflow and thus the wind starts with a very low bulk velocity which 
is why the mass-loss rate decreases. It is also observed that the sonic point of 
the accretion- 'self- wind' system is, in general, located closer to the event horizon 
compared to the outflow sonic point for all values of £, Msdd, lin and 70. 
Combining the informations obtained from the Fig. 2a and 2b, one can observe 
that for a flxed set of values of ^S,MEddilinilo^i the inflow sonic point (shown 
by Pg in the figures) and the outflow sonic points (for various 70) for flows in $2, 
is greater than those quantities obtained for flow in $1. Whereas the shock Mach 
number Mgh for $1 has a higher value compared to M^/i for <1?2- This is obvious 
because the potential $2 is steeper compared to $1 (see DS) hence $2 produces 
relatively higher acceleration and the accreting material gains more kinetic energy 
relatively faster to become supersonic at a distance relatively greater compared to 
the flow in $1 and as the sonic point of the inflow resides relatively further away 
from the event horizon, accreting material becomes supersonic at an earlier stage 
and the rate of increment of the ratio of dynamical (mechanical) to the thermal 
energy content of the accretion increases; also another contributing factor is the 
location of the shock, Eq. (15) also tells that for $2, in reality, should be less 
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compared to its value for $i. However, for any set of <^£, MEdd,7in,7oj , the mass 
outflow rate Rm for $i is less than that obtained in $2- We will come to this point 
in detail in next sub-section. 

3.3. Dependence of Rm on the specific energy of the flow 

In figure 3a, we have plotted the variation of Rm with the specific energy of the flow 

£ for a flxcd values of 7i„ (= |) and 70 (= 1.3) and for a number of values (shown 
in the figure) of the accretion rate (scaled in units of Eddington rate and shown in 
the figure as MEdd) of the inflowing material for flows in $1 (dotted lines), as well 
as for flows in $2 (solid lines). Any dotted line accompanying a solid line marked 
by a specific value of MEdd corresponds to the same accretion rate as that of the 
solid lines. We observe that for a fixed value of MEdd, lin and 70, the mass-outfiow 
rate non-linearly correlates with £ reason for which might be as follows: 
As £ increases, r^h decreases (see Fig. 1, §3.1) and the post-shock bulk velocity 
of the fiow u^y^ as well as the post-shock density increases. The outflow rate, 
which is the product of three quantities Vsh, pj'^ and zt^^ (see Eq. (lib)), increases 
in general due to the combined 'tug of war' of these three quantities. Moreover, 
the closer the shock is to the black hole, the greater the amount of gravitational 
potential will be available to be put onto the relativistic protons to provide stronger 
outward pressure and the closer the shock forms to the black hole, the higher is the 
post-shock proton temperature (the effective characteristic outflow temperature) 
and the higher is the amount of outflow (as the wind is observed to be strongly 
thermally driven, see discussion below). Thus the mass-outflow rate increases with 
£ because for a particular set of fixed values of MEdd, Rm is proportional to Mout 
which increases with £. 

The unequal gaps between the curves marked with different MEdd in the figure imply 
that when the inflow specific energy is kept constant, Rm non-linearly increases with 
the accretion rate of the infalling material. This is because, as £ is kept constant 
while MEdd is varied, the amount of infalling energy converted to produce the high- 
energy protons is also fixed. So the higher the value of MEdd, the larger the distance 
of the shock surface from the event horizon and the outflowing matter feels a low 
inward gravitational pull, the result of which is the non-linear correlation of Rm 
with MEdd- In other words, for a fixed value of 7i„, high energy high luminosity 
accretion produces more outflow for both of the pseudo-potentials considered here. 
Also one can observe that for a fixed set of parameters ^£,MEdd,7in,7o^, the 
amount of mass-outflow for <I>2 is higher compared to that in $1. 
In figure 3b, we plot Rm as a function of T^"^ for $1 (dotted line) and $2 (solid 
line) corresponding to the parameters used in the figure 3a. We see that post-shock 
temperature non-linearly correlates with the energy of the flow, and for a fixed 
accretion rate and adiabatic indices of the inflow and outflow, Rm also correlates 
with post-shock temperature, which indicates that the outflow is thermally driven 
as well. 
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3.4. Variation of Rm with adiabatic indices of the flow 

In previous cases, the polytropic index 7i„ of the accreting matter was always kept 
fixed at the value |. To have a better insight of the behaviour of the outflow, we 
plot Rm as a function of 7i„ (curves marked by 7i„ in Fig. 4) for a fixed value of £ 
(= 0.00001) and Msdd (= 1-0). The upper range of 7,„ shown here is the range for 
which shock forms for the specified value of £ and MEdd- We have chosen the value 
of 7o in such a way so that the condition 70 < 7j„ is always satisfied. Defining A-y 
to be A-y = 7j„ — 7o, we study the variation of Rm with 7j„ for A^ = 0.01. One 
can also obtain results for other values of A^ in the same way. Wc observe that 
Rm correlates with 7i„, which is expected bccaiisc the specific enthalpy of the flow 
increases with 7j„ to produce a higher post-shock temperature for higher value of 
7j„ (see also the curves marked by T^"^ in Fig. 4). We observe from our calculation 
that p^f^ and correlates while Vgh anticorrelates with 7;„. So increment of 7i„ 
satisfies all possible conditions to have a high value of Rm- In the same figure, curves 
marked by T^"^ show the variation of Rm with T^"^ corresponding to the values of 7j„ 
shown in the figure (scaled T^"^ — > x 10^-^^ to fit in the same figure) to show 
that here also the flow is thermally driven. Also, we observe that as well as the 
shock compression ratio Rcomp correlates with 7^ so that 'strong-shock' solutions 
are preferred to obtain a high value of mass-loss for this case. So we conclude that 
as the accretion approaches from its ultra-relativistic nature to its non-relativistic 
regime, the mass-loss rate increases. 

It is observed that for same set of parameters used to study the dependence of Rm 
on 7j„, flow in $2 produces more mass-loss compared to flow in $1, i.e., 



R 

for any value of 7i„ and the quantity 
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decreases with decrease of 7i„. However, in the figure, the curves showing the de- 
pendence of Rj^ with Jin are not properly resolved to explicitly show the difference 
for flows in $1 and $2 due to the scaling used here to show both the temperature 

dependence of Rj^ as well as the dependence of Rj^ on 7j„ in the same figure. 
We also study the variation of Rj^ with 70 (figure not presented in this paper) for 

fixed values of M^^d, 7m|- The general conclusion is that Rm correlates with 
7o for both of the pseudo-potentials. This is because as 70 increases, shock location 
and post-shock density of matter does not change (as 70 does not have any role in 
shock formation or in determining the Rcomp) but the sonic point of the outflow is 
pushed inward, hence the velocity with which outfiow leaves the shock surface goes 
up, resulting the increment in Rm- Here also we found that Rm for flow in $2 is 
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always greater than than Rm for $i and the deviation becomes more prominent as 
the value 70 approaches to 7i„. goes up, resulting the increment in Rm- 

4. Conclusion 

In this paper, we could successfully construct a self-consistent spherically-symmetric, 
polytropic, transonic, non-magnetized inflow-outflow system by simultaneously solv- 
ing the set of hydrodynamic equations governing the accretion and wind around a 
Schwarzschild black hole using two different pseudo-potentials proposed by Arte- 
mova et. al.^ Introducing a steady, standing, hadronic-pressure supported spherical 
shock surface around the black hole as the effective physical atmosphere which may 
be responsible for generation of accretion-powered spherical wind, we calculate the 
mass-outflow rate Rm in terms of only three accretion parameters (conserved energy 
of the flow £, accretion rate MecIcI scaled in units of Eddington rate and polytropic 
index of the flow 7i„) and only one outflow parameter (the polytropic index of the 
outflow, 7o). Not only do we provide a sufficiently plausible estimation of Rm, we 
could also successfully study the dependence and variation of this rate on various 
physical parameters governing the flow. 
The basic conclusions of this paper may be summarized as: 

1. Shock formation is not a generic phenomena, i.e., not all solutions contain 
shock, rather a specific region of parameter space spanned by £, Msdii and 
7i„ allows shock formation. For given values of £,MEdd and 7i„, while the 
value of shock location (in units of rg) correlates with Mscid, it anti- correlates 
with both £ and 7i„. High energy high luminosity purely non-relativistic 
accretion is a proper choice to produce high post-shock proton temperature ( 
thus to produce the maximum outflow). For other accretion parameters being 
the same, $1 produces a shock closer to the black hole compared to $2- This 
deviation is more prominent for high £ and 7i„. 

2. The shock surface can serve as the 'effective' physical barrier around the black 
hole regarding generation of mass loss via transonic spherical wind. The 
fraction of accreting material being blown as wind (which is denoted as Rm) 
could be computed in terms of three accretion parameters and one outflow 
parameter. 

3. Rm correlates with £,MEdd,lin and 70, Outflow could be generated for both 
sub-Eddington as well as super-Eddington accretion. For a fixed set of 

l^, MEdd,Jin,Jo^, mass outflow rate Rm is higher for flows in $2 compared 
to flows in <I>i, which indicates that the barionic content of the spherical wind 
is inversely proportional to the spatial gradient on the pseudo-potential used 
to study the problem though the exact physical reason behind this is not quite 
clear. 

4. If a shock forms, then whatever the initial flow conditions and whatever the 
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nature of dependence of Rm on any of the accretion/ shock/ outflow pa- 
rameters, Rrh normally correlates with post-shock flow temperature, which 
indicates that outflow is strongly thermally driven; hotter flow always pro- 
duces more winds. 

At this point, it is worth mentioning that the hot and dense shock surface around 
black holes, which is proposed here as the effective physical barrier around com- 
pact objects regarding the mass outflow, may be generated due to other physical 
effects as well for spherical accretion. Another very important approach 
launched recently was to construct such an 'effective barrier' for non-spherical disc 
accretion to introduce the concept of CENtrifugal pressure supported BOundary 
Layers (CENBOL). Treating the CENBOL as the effective atmosphere of the ro- 
tating flows around compact objects (which forms as a result of standing Rankine- 
Hugoniot shock or due to the maximization of polytropic pressure of accreting 
material in absence of shock), detailed computation of the mass-outflow rate from 
the advective accretion disks has been done, and dependence of this rate on various 
accretion and shock parameters has been quantitatively studied by constructing a 
self-consistent disk-outflow system. -^^'^^ 

Our calculations in this paper, being simply founded, do not explicitly include var- 
ious radiation losses and cooling processes, combined effects of which may reduce 
the post-shock proton temperature (which means the reduction of outflow tempera- 
ture), in reality could be lower than what we have obtained here and the amount of 
outflow would be less than what is obtained in our calculation. This deviation will 
be more important for systems with high accretion rates. Nevertheless, cases of low 
accretion rates discussed here would not be affected that much and our preliminary 
investigation shows that even if we incorporate various losses, the overall profile 
of the various curves showing the dependence of Rm on different infiow parameters 
would be exactly the same, only the numerical value of Rm in some cases (especially 
for high accretion) might decrease. 

Although in this work we have performed our calculation for a IQMq Schwarzschild 

black hole, general flow characteristics will be unchanged for black hole of any mass 
except the fact that the region of parameter space responsible for shock formation 
will be shifted and the value of Rm will explicitly depend on the mass of the black 
hole. 
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Figure Captions 



Fig. 1: Variation of shock location Vgh (plotted in logarithmic unit along Y axis) 
with the specific energy of the flow £ (plotted along X axis) for three values of 
accretion rate (scaled in the unit of Eddington rate MEdd) 0.25, 1.0 and 1.75 re- 
spectively from bottom to top. While the dependence for flows in $i is shown by 
dotted lines, flows in $2 is shown by solid lines; see text for details. 

Fig. 2: Combined solution topologies for transonic accretion-outflow system for 
flows in $1 (2a) and $2 (2b) for fixed values (shown in the figure) of accretion and 
outflow parameters. While Mach number is plotted along the Y axis, the distance 
(in units of r^) from the event horizon of the accreting black hole is plotted along the 
X axis in logarithmic scale. Solid curve marked with AB represents the pre-shock 
transonic accretion while four dotted curves represent the accretion-powered out- 
flow branches. Solid vertical line BC (marked by Vsh) with double arrow stands for 
the shock transition and solid line marked by CD stands for the 'self-wind' branch. 
Pg is the location of the inflow sonic point. See text for details. 

Fig. 3a: Variation of Rm with inflow specific energy £ for different values of 
accretion rate MEdd shown in the figure. Dotted curves correspond to fiow in $1 
and solid curves correspond to flow in $2- 

Fig. 3b: Variation of Rm with post-shock proton temperature T^"^ correspond- 
ing to values of £ shown in Fig. 3a. 

Fig. 4: Variation of Rm with the polytropic index of accretion 7j„ (curves marked 
by Jin) and the corresponding post-shock temperature (curves marked by T^f^) for 
flows in $1 (dotted curve) and in $2 (solid curve). The Mass-outflow rate non- 
linearly correlates with 7i„ as well as with post-shock proton temperature which 
indicates that the outflow is thermally driven as well, see text for details. 
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Fig. 2a: 
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Fig. 2b: 
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Fig. 3a: 
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Fig. 3b: 
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